x = (10:100:12501)';
y = lognpdf(x,8,0.5);
z = lognpdf(x,8,1.3);
stry = '\sigma = 0.5';
strz = '\sigma = 1.3';
y1 = 700;
y2 = lognpdf(y1,8,0.5);

z1 = 1800;
z2 = lognpdf(z1,8,1.3);

figure
plot(x,y,x,z);
set(gca, 'XTick', [0 3000 6000 9000 12000]);

title({'Two lognormal income distributions when \mu=8'})
set(gca, 'YTick', []);
text(4800,0.00012,stry)
text(180,0.00025,strz)
line([2626,2626],[0,0.00037],'Color','k','LineStyle',':')
line([2914,2914],[0,0.00037],'Color','k','LineStyle',':')


